Fake symmetry transitions in lattice Dirac spectra 
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In a recent lattice investigation of Ginsparg- Wilson-type Dirac operators in the Schwinger model, 
it was found that the symmetry class of the random matrix theory describing the small Dirac 
eigenvalues appeared to change from the unitary to the symplectic case as a function of lattice size 
and coupling constant. We present a natural explanation for this observation in the framework of a 
random matrix model, showing that the apparent change is caused by the onset of chiral symmetry 
restoration in a finite volume. A transition from unitary to symplectic symmetry does not occur. 

PACS numbers: 11.15.Ha, 12.38. Gc, 05.45.Pq 



I. INTRODUCTION 



In the past two years, Dirac operators T> satisfying the 
Ginsparg- Wilson (GW) condition Q, 



Z>75 + IsP = 2Vj 5 RD , 



(1) 



have attracted a great deal of attention in the lattice 
community because of their vastly improved chiral and 
topological properties. In a recent lattice study of the 
Schwinger model, Farchioni et al. Q investigated two dif- 
ferent versions of T>, a fixed point Dirac operator 2?f p j!J 
and Neuberger's Dirac operator X>n c @] ■ For the purpose 
of the present study, these two operators have essentially 
identical properties, and we concentrate on T>^ e in the 
following. The subtracted lattice condensate is defined 
by " 



(iplp) eab = --7 (trl?~ 

V \ ' gauge 



(2) 



where V = T>(1 — RD) 1 and V is the physical volume. 
The operator T> shares some important features with the 
continuum Dirac operator in Euclidean space. It is anti- 
Hermitian and satisfies {25,75} = 0. Thus, its nonzero 
eigenvalues occur in pairs ±i\ n with A n real. The spec- 
tral density of iV is defined by p(X) — Q^ n (A — ^n))/V. 
The Banks-Casher relation M then reads 



E = - lim (V>V>)sub 



lim lim 7ro(A) , 

A^O V->oo 



(3) 



where it is important that the thermodynamic limit is 
taken first. 

Our main concern in the present work are the some- 
what puzzling observations made in Ref. Q where the 
spectral properties of T> were compared to predictions of 
chiral random matrix theory (RMT) Q . RMT is a simple 
model which yields exact analytical results for the spec- 
tral correlations of the Dirac operator on the scale of the 
mean level spacing. Because of the chiral structure of the 
problem, one has to distinguish two regions in the spec- 
trum, (i) the bulk and (ii) the eigenvalues in the vicinity 



of zero (the latter is called the microscopic region or the 
hard edge of the spectrum). The analytical predictions 
of RMT are different in these two regions. Furthermore, 
there are three different symmetry classes, the chiral or- 
thogonal (chOE), unitary (chUE), and symplectic (chSE) 
ensembles H. 

Let us briefly summarize those findings of Ref. (|] 
which arc of relevance to the present work. The 
Schwinger model (QED2) is in the symmetry class of the 
chUE. For large volumes V cx L 2 /(3 (L is the number of 
lattice sites in each dimension and (3 = l/(ea) 2 is the di- 
mensionless coupling), it was found that the bulk as well 
as the hard edge of the spectrum are nicely described by 
the chUE predictions. As V decreases, the bulk prop- 
erties are still given by the chUE whereas at the hard 
edge, the data suggest a transition to chSE behavior. In 
Ref. Q , the sectors of topological charge v = and v = 1 
were investigated, and the apparent transition was seen 
in both sectors. These observations are puzzling, since it 
is not clear where the symplectic symmetry should come 
from. Also, the spectrum is not doubly degenerate as 
expected for the chSE. The authors of Ref. 0] already 
suggested that their observation might be an artifact of 
the small physical volume and that the agreement with 
the chSE may be accidental. 

In the framework of RMT, a transition between the two 
symmetry classes can be described by the simple model 



(1 - a) 2? c hUE + a £> c hSE 



(4) 



with Dirac operators of the appropriate symmetries and 
< a < 1 interpolating between the two ensembles. 
However, in such a model the transition from chUE to 
chSE is expected to occur simultaneously in the bulk and 
in the microscopic domain. We have confirmed this ex- 
pectation numerically. Thus, the simple ansatz (^) can- 
not explain the findings of Ref. Q . This is as it should 
be, since we do not expect a symplectic symmetry to be 
present in the first place. 

For RMT to be applicable at the hard edge, it is neces- 
sary that chiral symmetry is spontaneously broken, i.e., 
we require p(0) > 0. (Since on the lattice one is always 
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working at finite volume, we mean here the value ob- 
tained by extrapolating p(X) many level spacings away 
from zero to A = 0.) As we shall see, the apparent tran- 
sition from chUE to chSE symmetry in the microscopic 
domain is caused by the fact that p(0) vanishes as the 
physical volume decreases. This will be shown in more 
detail below, using a random matrix model which we con- 
struct in the following section. The spectral properties 
we consider are introduced in Sec. Ill, our results are 



presented and discussed in Sec. [V, and conclusions are 
drawn in Sec. |y| 



II. THE RANDOM MATRIX MODEL 



There are two parameters in the random matrix model, 
the linear lattice size L and the width 28 of the distri- 
bution P(ip) = 0(5 + ip)9(S — ip)/25. The parameter <5 is 
restricted to the interval < 5 < tt. It models the cou- 
pling [3 used in the lattice gauge simulation. Small values 
of 6 correspond to phases exp(iip) which fluctuate only 
weakly around unity. Thus, small values of 5 correspond 
to large values of f3 and vice versa. 

An analytical treatment of the problem is very dif- 
ficult. Therefore, we performed a numerical study for 
several combinations of the parameters L and 5. As we 
shall see, the random matrix model is able to reproduce 
and to explain the findings of Ref. Q] . 



Our approach closely parallels the construction of the 
Neuberger operator 0J. The Wilson Dirac operator for 
the d — 2 dimensional Schwinger model defined on an 
L x L lattice reads 

2 

M x . y = 5 x , y - \(1 - a /J ,)U tl (x)A Xty _ ll + 

(5) 

with hopping parameter k = 1/ (2m + 4) and Pauli ma- 
trices <7,. As in Ref. H, we will restrict ourselves to the 
case m = —1, k = 1/2. The boundary conditions for 
the fermions are periodic in space and anti-periodic in 
Euclidean time. This is taken into account by the func- 
tion A x>y ±f l which is the usual Kronecker delta 5 XtV ±f l 
except for the links from xi — L to X2 = 1 for which 
there is an extra factor of — 1. The U(l) gauge fields 
are represented by 2L 2 phases U^(x) which obey periodic 
boundary conditions and fluctuate in the update process. 
In our model, they are simply replaced by independent 
uncorrelated random phases exp(iy>) with tp drawn at 
random from the interval [—6, 5]. For a single flavor with 
two spinor indices, M. thus becomes a (sparse) random 
matrix of dimension 2L 2 . 

From the matrix W obtained by replacing the gauge 
fields by random phases we construct the Neuberger op- 
erator 01, 



W Nc = 1+75 e(7sW) 



(6) 



where e is the sign function e(A) = AjvAM. 2?Ne 
(and hence also Wnc) satisfies the GW condition (jl]) with 
R = 1/2. Its spectrum is located on the unit circle in the 
complex plane with center at z = 1. Furthermore, for 
certain background fields it possesses exact zero modes 
whose number is related to the topological charge of the 
background field. In analogy with the definition of T> we 
finally construct the random matrix 



W = W Nc (l->VNe/2)- 



(7) 



This projects the spectrum of W^e onto the imaginary 
axis. In the following, we will investigate the spectrum 
of the matrix W. 



III. SPECTRAL PROPERTIES 



Since our main purpose in this work is to understand 
the puzzling observations of Ref. ||, we have investi- 
gated the same quantities considered therein. These are 
the spectral density p(X), the distribution of the smallest 
(positive) eigenvalue P(A m ; n ), and the nearest neighbor 
spacing distribution P(s). The quantity s is the spacing 
between adjacent eigenvalues in the bulk after the spec- 
trum has been unfolded ||. P(A m i n ) and P(s) measure 
the spectral properties at the hard edge and in the bulk 
of the spectrum, respectively. (We will not discuss the 
microscopic spectral density which was also considered 
in Ref. 0] since for the purpose of the present study, 
-P(A m in) is sufficient to reveal the properties of the spec- 
trum near A = 0.) The global density p(A) of the random 
matrix model is not a universal quantity, i.e., it depends 
on the details of the model and is not expected to agree 
with real lattice data. However, -P(A m ; n ) and P(s) are 
universal for a specific symmetry class, i.e., these func- 
tions are insensitive to the details of the dynamics, and 
the results obtained in RMT should agree with those of 
lattice simulations (unless the conditions are such that 
RMT is not applicable, see below). 

The nearest neighbor spacing distribution P(s) in the 
bulk of the spectrum is expected to agree with the RMT 
result irrespective of whether or not chiral symmetry is 
broken In other words, it is insensitive to the value 
of p(0). P(s) can be unambiguously constructed from 
the lattice data without any free parameter by unfold- 
ing. The RMT results for P(s) are well approximated by 
the Wigner surmise B 



poo 



(32/7r 2 )s 2 exp(-4s 2 /7r) chUE 
(2 18 /3 6 tt 3 ) s 4 exp (-64s 2 /9tt) chSE . 



(8) 



The nearest neighbor spacing s is expressed in units of 
the mean level spacing in the bulk. Note that P(s) does 
not depend on v or on the number of massless flavors Nf. 

In contrast to P(s), P(A m i n ) is only given by the RMT 
result if chiral symmetry is broken, i.e., if p(Q) ^ 0. Also 
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in contrast to P(s), it depends on the number of zero 
modes v and on Nr. In the following, we restrict our- 
selves to the quenched case, Nf = 0. The energy scale 
for a comparison between the RMT result and numer- 
ical data is set by the mean level spacing at the hard 
edge, 1/Vp(0) — n/VT,, which can be determined by 
Eq. (||) without resorting to RMT. In terms of the vari- 
able z — FSAmin, the RMT results for the chUE are 
given by [jlO) 

f (z/2)e-* 2 / 4 i/ = 0, 

P(z)= \ (z/2)e- 2 / 4 / 2 (z) u = l, (9) 

[ (z/2)e~* 2 ^[ll(z)-h(z)I 3 ( Z )] y = 2, 

where I denotes the modified Bessel function. The RMT 
results for the chSE are more complicated and can be 
found in Ref. plj . The only simple case is v = for 
which 111 



P{z) 



'■e- z2 ^I 3/2 (z) 



(10) 



As we shall show in Sec. |v|, the apparent change from 
chUE to chSE symmetry in P(A min ) is caused by the fact 
that p(0) (i.e., the extrapolated value) vanishes in a finite 
volume for sufficiently weak coupling. An unambiguous 
and parameter-free comparison of numerical data with 
RMT results is only possible if p(Q) and hence E are 
nonzero. In the present work, however, we are particu- 
larly concerned with the region in the parameter space 
(L, 8) where S — > 0. In this case, we are faced with the 
problem of how to set the energy scale for A min . We cir- 
cumvent this problem by employing three independent 
methods of analyzing the data. 

(1) We rescale the smallest eigenvalue, both for the 
numerical data and for the RMT results, such that 
(Amin) = 1- In this way, the various cases we consider are 
characterized only by the shape of the distribution; the 
ambiguity in the determination of the scale is avoided. 

(2) For each data set, we construct the ratio 



V (A min ) - (Amin)" 



(11) 



which eliminates the energy scale. The RMT results for 
this ratio are 



(12) 



(3) We fit the data for P(A min ) to the RMT results 
of both the chUE and the chSE (for fixed v) and apply 
goodness-of-fit tests |uj. For the determination of the 
scale, we use the method of maximum likelihood. For 
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the quality of the fit, we use three tests: 
(3a) A Smirnov-Cramer-Von Mises test of the unbinned 
data (see Sec. 11.4.1 of Ref. This test yields a num- 

ber NW 2 which for a good fit is smaller than one (N is 
the number of configurations). 

(3b) A chi-square test of the binned data using equiprob- 
able bins, with an "optimal" number k of bins scaling like 
N 2 / 5 (see Sec. 11.2.3 of Ref. g§). For a good fit, the re- 
sulting x 2 /(k — 1) is not much larger than one. 
(3c) The probability Q that the chi-square should exceed 
the value of \ 2 computed in (3b) by chance fl4| . This 
is given by Q = 1 — P((k — l)/2,% 2 /2) with the incom- 
plete gamma function P(a, x) = f Q dtt a ~ 1 e~ t /T(a). If 
the probability Q is very small, the theory can be statis- 
tically rejected. 

Let us emphasize again that in the region where chi- 
ral symmetry is spontaneously broken, P(A mm ) is un- 
ambiguously given by RMT without any free parameter. 
Here, however, we explicitly consider the transition re- 
gion where p(0) 0. We use the somewhat more in- 
volved methods (l)-(3) to make clear and parameter- 
independent statements on the behavior of -P(A m in) in 
this transition region. 



IV. RESULTS AND DISCUSSION 

We have scanned the parameters L and 8 to reveal the 
main properties of the spectrum. For each parameter set 
we generated several thousand random matrices W as de- 
scribed in Sec. [il]. For each realization the full spectrum 
was computed, and the configurations were grouped ac- 
cording to their number v of zero modes. The quantities 
we consider were then averaged over realizations of W. 

In Fig. |l|, we have fixed the parameter L = 10 and var- 
ied the parameter 8. We have plotted the spectral den- 
sity p(A) in the vicinity of A = for all v combined, and 
the rescaled distribution of the smallest positive eigen- 
value -P(Amin) for v — 0, 1, 2. (Note that the zero modes 
are not included in p(X).) Also given in the figure are 
the values of the ratio r, see Eq. (|ll"|). The results of 
the goodness-of-fit tests (3a)- (3b) defined in Sec. Ill are 
summarized in Table ||. 

The main message of this paper can be read off im- 
mediately from these plots. The largest value of 8 = 1.6 
corresponds to the strongest coupling and hence to the 
largest physical volume. The extrapolated value of p(0) is 
clearly nonzero. Thus, RMT is applicable, and P(A m in) 
agrees very well with the chUE curves for all three val- 
ues of v. The data also pass the other tests proposed in 



Sec. Ill: The ratio r agrees well with that of the chUE, 
see Eq. (|l2|), and the goodness-of-fit tests, see Table §, 
show that the chUE describes the data well. 

On the other hand, the smallest value of 8 = 1.2 corre- 
sponds to the weakest coupling and hence to the smallest 
physical volume. In this case, the extrapolated value of 
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FIG. 1. Global spectral density p(X) (left column) and rescaled distribution of the smallest eigenvalue P(A m i n ) of W for 
L — 10 as a function of S for v — 0,1,2. The solid and dotted curves in the plots for P(A m i n ) are the chUE and chSE 
predictions, respectively. The quantity r is defined in Eq. dill). 



p(0) vanishes, and chiral symmetry is restored. Thus, a 
comparison with RMT becomes meaningless in the mi- 
croscopic region. The shape of the distribution P(A m i n ) 
as well as the ratio r of Eq. (|llj) are now different from 
the RMT predictions. Neither the chUE nor the chSE 
provide good fits to the data, as can be seen from the 
goodness-of-fit tests in Table |. For the intermediate case 
of S = 1.3, the value of p(0) is neither clearly zero nor 
clearly nonzero. Figure |l| and Table | indicate that al- 
ready in this case RMT no longer provides a good de- 
scription of the data. 

Equivalent conclusions can be drawn from Fig. |^ in 
which we exhibit similar plots, but now for a fixed value of 
5 = 7r (corresponding to "strongest coupling" ) and L var- 
ied from 10 to 4. The goodness-of-fit tests are also shown 
in Table |. Again, as the physical volume decreases, the 
extrapolated value of p(0) vanishes, and the agreement 
between numerical data and RMT breaks down. We also 
observe from Fig. ^ that RMT ceases to describe the data 
first for the nonzero values of v and then for v = 0. 

A remarkable fact, however, is that there are some 
combinations of the parameters L and 8 for which the 
shape of the distribution P(A m i n ) resembles the chSE pre- 
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TABLE I. Results of the goodness-of-fit tests as explained 
at the end of Sec. [II, Zero means zero to machine precision. 



4 



0.15 
p(X) 
0.1 




J 1 I I 1 1 , I 1 1 , 


1 1 1 1 1 '_ 


- L=4 








Hi.. 


, , , 1 , , " 



- r = 3.60 


i/=2 - 


xJ<y\ 





0.5 1 , 1.5 2 

A . 




1 1 1 1 1 1 1 1 

- r = 5.68 1 

L I 


i i i i i i i i_ 







0.5 1 , 1.5 2 0.5 1 . 1.5 2 



FIG. 2. Same as Fig. hi but keeping S = n fixed and varying L. 



diction. This effect is particularly pronounced for the 
nonzero values of v, which was also found in Ref. j^] (for 
v = 1). In our data, this is seen in Fig. |l| for S = 1.3, 
v = 1,2 and 5 = 1.2, v = 2, respectively, and in Fig. |2] for 
L = 4, v = 1. (Recall that the plots of -P(A m i n ) are not 
fits but are obtained by rescaling such that (A m i n ) = 1.) 
However, neither the ratio r nor the goodness-of-fit tests 
in Table i] indicate that the chSE really provides a quan- 
titative description of the data. The fact that the chSE 
seems to work better than the chUE should not be taken 
too seriously since neither of the two ensembles is sup- 
posed to be applicable in these cases. The apparent 
agreement of the shape of the curve with the chSE is, 
for all we know, purely accidental and should not be in- 
terpreted as a transition to symplectic symmetry. 

We have also constructed P(s) for the same parameter 
values that were used in Figs, [j] and [|. Two examples 
are shown in Fig. |3|. In agreement with the observations 
of Ref. [||, we find that P(s) is always described by the 
chUE prediction, even in those cases where p(0) vanishes. 
This is consistent with the fact that, even in the chirally 



symmetric phase, the bulk spectral correlations of the 
Dirac operator are still described by RMT M . 
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FIG. 3. Nearest neighbor spacing distribution P(s) for 
some of the parameter values for which P(A m i n ) is not de- 
scribed by the chUE. The histogram represents the data, and 
the solid and dotted curves are the chUE and chSE predic- 
tions of Eq. M , respectively. 
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V. CONCLUSIONS 
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Let us summarize the picture which has emerged from 
the study of the random matrix model. The Dirac oper- 
ator of the Schwinger model has the symmetries of the 
chUE, and our model has the same symmetries, the only 
essential difference being the replacement of the U(l) 
gauge fields by random phases. The bulk spectral cor- 
relations of the Dirac operator, measured by P(s), are 
described by the chUE result for all values of the param- 
eters, regardless of whether or not p(0) is nonzero. This 
is consistent with the findings of Ref. The micro- 
scopic spectral correlations, here measured by P(A m ;„), 
are given by the chUE predictions only for those val- 
ues of the parameters for which p(0) > 0. By changing 
the parameters one can decrease the physical volume and 
cause p(Q) to vanish. The small eigenvalues are then no 
longer described by RMT. In this case, however, there 
exist some special values of the parameters for which the 
shape of the distribution of the smallest eigenvalue re- 
sembles that of the chSE. However, other measures such 
as the ratio r of Eq. ( |il| ) and the various goodness-of- 
fit tests clearly show that the chSE does not describe 
the data. Thus, the apparent agreement with the sym- 
plectic symmetry case is an artifact of chiral symmetry 
restoration in a finite volume and should be regarded as 
accidental. 

There is a simple lesson to be learned from our analysis. 
If one wants to use random matrix methods to analyze 
data from lattice simulations, one has to make sure that 
one is working in a regime in which RMT is applicable. 
For the bulk spectral correlations, this is never really an 
issue. However, for RMT to describe the small Dirac 
eigenvalues it is necessary that p(0) > 0, i.e., that chiral 
symmetry is spontaneously broken. 
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